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Abstract 

Exact massive S-matrices for two dimensional sigma models on symmetric spaces SU (2 N) /Sp(N) 
and Sp(2P)/Sp(P) x Sp(P) are conjectured. They are checked by comparison of perturbative 
and non perturbative TBA calculations of free energy in a strong external field. We find the 
mass spectrum of the models and calculate their exact mass gap. 


1 Introduction 

Recently renewed interest to two dimensional sigma models with and without topological theta 
terms on different coset and, in particular, symmetric and supersymmetric spaces was stimulated 
by their relevance for Integer Quantum Hall effect [1] - [6], and two dimensional fermionic systems 
with different types of symmetry and quenched disorder (see, for example,[7] - [10],[2]). This is 
the reason why any new exact solution of such models, being important by itself in general, can 
also find a practical application. 

Essential part of symmetric space sigma models are exactly solvable, since were proved (or 
believed) to be integrable. Classical integrability of sigma models on any symmetric spaces G/H 
is known for a long time [11], but quantum integrability for symmetric space with H non simple 
may be destroyed by possible anomalies [12]. But even for H simple, when quantum integrability 
is undoubted, exact S-matrices are known not for all symmetric spaces. (For Cartan classification 
of symmetric spaces see for example [17]). 

Recently [13] the list of known exact S-matrices for integrable sigma models on symmetric 
spaces was extended by new massive and massless S-matrices for spaces AI ( SU(N)/SO(N) ) 
and BDI ( 0(2P)/0(P) x 0(P ) ) without and with theta term. The conjectured S-matrices 
were checked by comparison of T = 0 TBA calculations of free energy of the system in strong 
external field with perturbative calculations [13], and also by T —> oo TBA extraction of UV 
central charge for these two models [14]. It is remarkable that BDI sigma model turns out to be 
quantum integrable in spite of possibility of anomalies. 

There are other close relatives of the models considered by Fendley: sigma model on All ( 
SU(2N) / Sp(N ) ) is close to AI, and sigma model on CII ( Sp(2P) / Sp(P) x Sp(P) ) is analogous 
to BDI sigma model - in both cases orthogonal (sub)group is changed to symplectic. We are 
going to show that the analogy extends also to S-matrix conjecture, providing one more example 
(CII) of quantum integrability of non supersymmetric sigma model on a factor group manifold 
with non simple invariant subgroup. 

*address:Shahal str. 73/4, Jerusalem, 93721, Israel; e-mail:ababichenko@hotmail.com 


1 



The plan of the discussion is the following. The first section describes standard perturbation 
theory technique of free energy calculation in a strong external field. In the second section 
we conjecture our fundamental S-matrices for the models using some symmetry arguments in 
support of it, and show how the mass spectrum of the model follows from bootstrap. After 
that we calculate the same free energy using T —► 0 TBA technique based on the conjectured 
exact S-matrix. We show correspondence between perturbative and non perturbative results 
which confirms correctness of the conjectured S-matrices. Moreover, two expressions for the 
same quantity - the free energy - fixes the mass gap for the models exactly. We conclude by 
brief discussion of results and give some technical details for perturbative calculations and the 
S-matrices, like explicit form of projectors, in the Appendix. 

2 Perturbative analysis 


We are considering the Lagrangian 

C 0 = ^Tr{d^d^} 

where we introduced a matrix representation d> for coset G/H elements. Important property of 
this Lagrangian is its global G invariance. We work with Lie algebraic current representation of 
Lagrangian: 

C 0 = -^Tr{(g~ 1 d tl g)(g- 1 d IJ ,g)} (1) 

and consider representation of a coset group element as an exponentialization of its Lie algebra: 
g = exp(i J2i niE 1 + i [C-, where H-j are generators of Cartan subalgebra of a coset space, 

and E 1 - other its generators. Explicit form of the basis Ej we use for calculations one can find 
in the Appendix. 

In order to check our conjectures about an S-matrix of the model (see the next section), 
we do, for today, a standard procedure of putting the system into a strong external field [15]. 
Strong field (and hence high energy limit) gives an opportunity to believe to perturbation theory, 
because of asymptotical freedom of the models. So we replace derivatives in (1) by ’’covariant” 
ones: 

D,,g = d tl g - hS lt0 (Qg + gQ) 

where h is a strength of the external field, which is chosen in the direction if in the Cartan 
subalgebra of G: Q = ( hMi = (qH). The lagrangian density (1) in the presence of the 

source becomes 

C = C 0 + ^ Tr{(g~ l Q + Qg~ x ) d 0 g} - Tr{Q 2 + g^QgQ} (2) 

We are going to calculate the dependence of the free energy on the external field h: 5f(h) = 
f(h) —/(0) using perturbative calculation in the running coupling constant A (h). We will restrict 
ourselves by quadratic part of the euclidian lagrangian in the fields raj, which turns out to be 
enough for our purposes. Some details of calculations one can find in Appendix and the result is 


£~-^a 2 + -£' 

L, — ^ q i ^2 ^0 
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( 3 ) 


where for both All and CII cases 

1 N 

£, o = 9 ( n b((^) 2 + tfMijiq))™*] + mijddf,) 2 + /i 2 M^(g))m- i } 

i>3 =i 

with mass matrices 

M[j(q) = (<fe — <Zj) 2 + (qi+N ~ Qj+N) 2 
M-f(q) = (■ qi - q j+ Nf + (qi+N ~ qjf 

for All case and 

M ij(q) = (qt - qj+p ) 2 , M ij (q) = (% + qj+p ) 2 

for CII case. (Here in (3) there is no need in complex conjugation * in CII case.) 

Lets point out that at this level some of the fields nj decoupled, and we wrote only those which 
enter in an h dependent manner. At the tree level one has Sf(h )o = — p-(/ 2 , q 2 = Y%=i 2P ^ <7?■ 
Free energy at the one loop level is just properly regularized \ l n det((9 /i ) 2 + h 2 Mij (q)). Using 
standard dimensional regularization e = d — 2 one finds 

6f{h)l = ~l^F q2 + E E^t^) 1-7 E + ln4vr - ln (^M^q)^ 

where fi is a mass parameter of dimensional regularization, and Piq 2 = J2i J2i>j =l Afl (g). After 
some algebra one can find that P\ = 2N for All case, and (5\ = 2P + 1 for CII case. Here we 
use the condition ^iV( 2 P) ^ necessary in All case, and just will correspond to our 

concrete choice of the external field in CII case (see below). 

The point is that the quantity Sf(h ) is renormalization group invariant when A runs with /i, 
so we can set fi = h and use the results of /3-function calculations (without external field), done 
for almost all symmetric spaces [21] up to three loops. We need the result up to two loops: 

^(A 2 ) = -i-(AA 4 + /3 2 A 6 ) - 0(A 8 ) (4) 

where /3 2 = 2N(N — 1) for All and /3 2 = P(2P + 1) for CII , and Pi is the same as above, 
since our calculation reproduced the correct form of one loop beta function. So adding necessary 
counterterm to lagrangian we get the following expression for the free energy 

Ah 2 h 2 

$f(h) = -jnxiV 2 - J2 M ij(q)i^ M !j(q) -1] + 0(A 2 ) (5) 

' I i>j= 1 

One can solve equation (4) 

1 „ . h /3 2 h ( 1 

—p— = Pi ln -1- p- In ln -b O - j— 

\ 2 (h) A ms Pi A-ms \lnpp 

where A ms is the cutoff parameter of minimal subtraction scheme, and substitute it into (5): 

sm = - 4 ' ,2 " 2A ( ln + | 1ln 1,1 + c ) (6) 

c = E E (7) 

J i>j = 1 

This expression will be used for comparison with result of free energy calculation by TBA based 
on exact S-matrix, which we are going to present now. 
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3 Exact S-matrices 


As in many other examples of quantum integrable models with higher rank Lie algebraic (ac¬ 
tually Yangian) symmetries, one can expect that particles group into multiplets corresponding 
to irreducible representations of the symmetry. As we mentioned, there is a global G symmetry 
acting on the coset space G/H, so we assume the S-matrix is related to branching rules of decom¬ 
position of highest weight reps of G into irreducible reps of H. The fundamental S-matrix usually 
is related to the shortest highest weight reps of G. In [1] a general matrix form was suggested, 
to which one can transform any factor group element of symmetric spaces by choice of a proper 
H gauge. These forms are quite useless for us here, but one can see that these matrix forms may 
be done antisymmetric for both All and CII cases. The minimal antisymmetric representations 
in su and sp algebras are reps with highest fundamental weight / 12 . This gives rise to our con¬ 
jecture: the fundamental S-matrix of All and CII symmetric space sigma models are described 
by rational x p , 2 S-matrices of SU ( 2 N) and Sp(2P) symmetry correspondingly. As it is well 
known, Lie algebraic symmetry with crossing and unitarity does not fixes S-matrix completely. 
The remaining so called CDD ambiguity is very important, since it in particular may change the 
pole structure of the S-matrix, i.e. defines bound states and spectrum of the model. This CDD 
ambiguity should be resolved using any kind of arguments, e.g. physically required coincidence 
of the S-matrix to a some known one, at a specific value of one of its parameters. There are 
two types of rational S-matrices of general series of Lie algebraic symmetries. Gross-Neveu like 
S-matrices have additional CDD factors with poles which, through the bootstrap, lead to a set 
of massive multiplets corresponding to all highest fundamental weights reps, while sigma model 
like S-matrices usually are ” minimal” (have no poles in the physical strip of rapidity) and have 
no these CDD factors. As it was conjectured and confirmed by different checks [13] [14], sigma 
models on the symmetric spaces AI and BDI , are similar rather to Gross Neveu models, since 
they have bound state coming from CDD poles. The same happens in our case. 

The fundamental S-matrix of All sigma model has the form 


S m (0) = X(x)s min (x ) (V 2/i2 + + t-x2A-x ^) (8) 

with 

sinh \ (6 + 47ri/21V) _ sin7r(2A + x) 
sinh \(9 — 4m/2N) sin7r(2A — x) ’ 
o M _ r(l - x) T(A + x) T(1 — A — x) T(2A + x) 
mml J T(1 + x) T(A — x) T(1 — A + x) T(2A — x) U 

where x = - 577 , A = and T^-projector on a rep with highest weight to. Explicit form of the 
projectors one can find in Appendix. 

For the CII sigma model we conjecture the following form of the fundamental S-matrix 


S m (x) = X(x)S min (x) (v 2 , 2 + \_ *. + _ X V 2lU1 + l _ ^ - 1 3 - V 0 (10) 

\ 2 / » X 2 • X 2 X 

A + A + x A — A + a; A + x 2A + a; \ 

+ 7T-+ A-T-X-+ A-TA- V W 

A — x A — - A — x A — x 2A — x I 

. _ A + x T(—A — x) T(2A + x) T(1 — x) T(A + x) 

x ) - A _ x r ^_ A + “^ r ^ 2A _ + r ( A _ x 
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X{x) 


w r(i + x ) r(i + a - x) r(i - a + ®) r(± + 2A - x ) 
x r(i - x) r(i + a + x ) r(i - a - x ) r(± + 2A + x ) 

sin tt(x + 2A) s i n71 ( x + 2 — ^a) 
sinvr(x - 2A) sin ^ (x - ± + 2A) 


where A = 2(2 ^ +1) . 

In both cases the product 5 m j n and parenthesis is the minimal unitary and crossing symmet¬ 


ric S-matrices of SU{2N) and Sp(2P) symmetry, fused from corresponding elementary vector 
representation S-matrices [18] [19] [20]. They have no poles on the physical strip of rapidities 9. 
Additional CDD factors in both cases provide the only pole (and hence a bound state particle) in 
the V IH channels (in Sp case there is also cross channel pole). It is clear that S-matrix describing 
the scattering of particle from fundamental P 2 multiplet on the particle from ^4 multiplet may be 
obtained by fusion and will give a pole in projector, and so on. In this way we get a spectrum 
for both models described by p, 2 k multiplets,( k = 1,..., N— 1 for All case and k = 1,..., P for CII 
case). Mass spectrum can be written as mk = M sin (^ 77 ) for All , and mu = M sin 
for CII, where M is a mass scale. 

As it sometimes happens in series of higher rank symmetric integrable models, in their lowest 
rank cases they often coincide with some other series (for example, lowest rank thermally per¬ 
turbed WD n CFT n = 2 are just parafermions with a proper perturbation). A remarkable hint 
for integrability of the models we are considering here, we get from the fact that SU {A) / Sp{2) 
is isomorphic to SO(6)/SO(5). It means that at N = 2 our All S-matrix should have the 
well known form of 0(6) sigma model. In the same way for CII case Sp{2)/Sp{l) x Sp{ 1) ~ 
SO(5)/SU(2) x 577(2) ~ 5'0(5)/50(4) - it is the 0(5) sigma model. Lets recall that 0(K ) sigma 
model has only one (vector) multiplet of particles in the spectrum (no bound states) with the 
following S-matrix 

oO(K),o\ T (1 — x)T{\ + x) T(x + 7 4 ^)T(i +- x) ( x +7^2.7-, 

5 w = UM , 77F7T-— 1 uvi —I —~ \ V s + Z - ~ Va 


+ 


± T \2 - i T Tc^2) l _r TTZ ^ 


K—2 


X — 


After some T function algebra one can see that K = 6 case really coincides with (8) at N = 2 , 
with corresponding representations mapping Vs —■► Vi^^A 'Pvi+w'P 0 “^ • Ti the same 
way one can check that K = 5 S-matrix is the same as (10) at P = 1. In this case representation 
correspondence has the form V s — > —*- Vim ■ Vo — > Vo ■ (Projectors V^+^V^V^ are 

absent in this case.). 


4 TBA calculation of the free energy 


The main point in T = 0 TBA analysis of our models in external field is based on a skill to chose 
external held in such a way, that TBA system will be the simplest, i.e. the ground state will 
contain the minimal number of particles generated from vacuum by external held. The fact that 
the held is strong gives a basis for the assumption that only particles with maximal charge will 
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be generated by external field. As we said the fundamental S-matrices for both All and CII 
models are rank two antisymmetric tensors and an external field from the Cartan subalgebra 
of G ,A = diag 2 N, 2 N{ A\. ...,A 2 n} ,acts on them as Aaij = (A % + Aj)aij. 

For the All case we chose the field in the form 

A = diag 2 N,2N{l, 1, ~ N _ ^ -> ~ N _ f ) l 11 ) 

(We work with the same normalization for fields and charges, and the meaning of our normaliza¬ 
tion choice -4= will be clear below). Then the ground state will contain only particles of the type 
ai 2 , since they have the maximal charge 2. One can see using the explicit form of projectors (see 
Appendix) that the scattering process ai 2 + <212 —► 012 + <212 takes place only in the P-2^ channel 
and the S-matrix for it is a prefactor before the parenthesis in (8). 

The situation is more complicated in CII case. We chose 

A = -y=diag4P } 4p{l, —1,0,..., 0, -1,1, 0, ...0} (12) 

with non zero elements on the places 1, 2, 2P +1, 2P + 2. In principle any combination of particles 
which is 0(2) invariant and has a maximal charge in the field A, can serve as a representative 
of the ground state. One can see that the combination d = a^p +2 — a2,p+i + ai,p+i — 02,P+2 
has the maximal charge 2. In addition this particle has two important properties, which one can 
find analyzing the projectors (see Appendix):firstly, (dd)ij = 0, and hence P /t2 . P 2 /tl ■ Pq are zero 
for the scattering of d on itself, and , secondly, this particle scattering on itself does not produce 
other particles and amplitude of the scattering is the coefficient before the projector P 2/i2 - the 
prefactor before the parenthesis in (10). 

So in both All and CII cases, with the choice of external field we described above, we have 
one particle in the ground state. Following standard technology of thermodynamic Bethe ansatz 
(TBA), one can get the TBA equation for the so called dressed energies of the particles 

r B 

e{6) = h — m cosh 6 + / dd'(f){9 — Q')e{Q') (13) 

J-B 

in terms of which free energy as a function of external field is 

m r B 

F{h)-F{ 0) =-/ d6 cosh{6)e{6) (14) 

27r J-b 


Here B is a function of h/m determined by the boundary condition e(±P) = 0, and 0 is a kernel 
defined by the S matrix S(9) = X(Q) S m m (9) for the scattering of the particles: 


^ = h^ s{0) 

The Wiener-Hopf method of solution of the integral equation of the form (13) in the limit h/m 
00 , aimed to extraction of the free energy (14), gives an answer for it 


K 


F(h)-F(0) = --h 2 


. h 1 h , 

In-1- (s + -) lnln- he + ... 

777 . 2 m 


(15) 


if Fourier transform of the kernel (f>{ul) = 1 — K(uj) factorizes into K(uj) = y, where 

K±(u>) are bounded and have no poles or zeros in the upper(lower) half plane and have an 
asymptotic for small £ 


K+m 


K 

7 ! 


e -^ in « 


(1 — 6 £ + 0 (£ 2 )) 


(16) 
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with some constants s, k and b, and 


/ i V^e b 

c = in--—--1 + 








(17) 


The detailed proof of this statement one can find in [22], Calculation of the Fourier transform of 
the kernels gives 


K(u) = 2e ,r|w|A 


sinh(7r|w| A) sinh(7ru;(l — 2A)) 
sinh(7ru;) 


(18) 


for All case and 


K(u)) = 


sinh(7r|w| A) cosh(7rcj(| — 2A)) 
cosh(7ro;/2) 


(19) 


for CII case.By |u;| we mean here the function which has the following analytical continuation 
to the whole complex plane: |cu| —> u * sign(Re(uj)). Factorization of (18) may be done as 


AAM = 


/ A(1 -2A)o; T(-iujA)T(-iuj(l - 2A)) 

27 t r(-iw) 


where the boundness of K + requires /i = A In A + (1 — 2A)ln(l — 2A), which leads to the 
asymptotic behavior of the type (16) with the constants 


s 


—A, K 


1 

i/27tA(1 - 2K) 


b 


M - A 'fE 


Kernel (19) factorizes as 


K + (u) 





r(—?'wA)r(i — 

■p/l_ ioj \ 

1 V2 2 / 


2A)) 


with /i = Ain A+(i —2A) ln(i —2A) + i ln2. It has asymptotic (16) with the following constants: 


s = —A, K = 


V2ttA 


, b = n + jeA 


5 Comparison of results 

We are going now to compare the results of perturbative (6) and non perturbative (15) TBA 
calculation of the free energy. First of all, using the form of external fields (and hence the 
charges) we chose (11),(12), one can see that prefactors of parenthesis coincide for both All and 
CII cases. Although the normalization of charges is ambiguous total factor, the fact that in this 
prefactor we get in both cases the same dependence on N and on P, is highly nontrivial check 
of not only our S-matrix conjecture, but also of the conjecture about the particle content of the 
ground state. Second, even more impressive, check is the comparison of the coefficients before 
the subleading term lnln h in both formulas. One of them is defined by beta function coefficients, 
another - by purely exact S-matrix dependent TBA analysis. Again, they coincide for both All 
and CII cases. Moreover, after we saw the coincidence of the leading and subleading terms in the 
limits h » A-ms, h » m, one can use different expressions (6) and (15) for the same quantity in 
order to fix the relation between mass scale m and the renormalization scheme parameter A ms- 
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This comparison involves the constant terms c, c' in both expressions. In the leading order of big 
h one just has what is called exact mass gap for the models 

777 

In -—— = c — d + 0(1/ In h ) 

A MS 

Using the form of the external field (11),(12) and the expressions for c, c/(7),(17), one can calculate 
the gap explicitly for both sigma models. 

6 Discussion 


We proposed fundamental S-matrices for two dimensional sigma models on All and CII sym¬ 
metric spaces. We checked them by comparison of perturbative and TBA calculations of free 
energy in a strong external field of a specific form and found the desired correspondence. In 
this check not all the particles of the conjectured spectrum have participated - only subsector 
of fundamental S-matrix was used. In this sense T > 0 TBA check based on extraction of UV 
central charge seems to be a more complete, since it involves all the particles of the spectrum. We 
hope to report on this soon [24], 

The quantum integrability of All sigma model was expected since the factorization subgroup 
is simple in this case, but the integrability of CII sigma model is a ’’surprise” in this sense, 
because one might expect anomalies for non local current conservation [12], It is clear that a 
deeper understanding of their integrability is desired from the point of view of conserved currents 
algebra and their symmetry. More exactly, the question remaining unclear is what kind of Yangian 
symmetry is responsible for integrability of sigma models on coset spaces. For today, the only 
known to us mathematically rigorous formulation of coset like symmetric Yangians are twisted 
Yangians, but they are known to be responsible for boundary integrability of sigma models [23], 
and hardly have something to do with quantum integrability without boundary. 
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8 Appendix 

8.1 Bilinear action 

If we chose the symplectic form as 2N by 2N matrix of the following block form ^ ^ 

basis (non orthonormal) for the generators of Lie algebra of All symmetric space one can chose 
in the form of the following 2N by 2N matrices 

— i E ij + E J+ N )t +N ), * 7 ~ j 

— (■A'ij+JV Ej,i- |-jv) 

— (A/j+ATj Ej-\-N,i) 

= (Eii — Ej+iy+i + E i+ N,i+N — A'j_|_7v_|_i i i_|_7v_|_i) 


E\j 

E iJ 

pin 

Eij 

Hi 



where i,j in the first three lines are running from 1 to N, and in the last (Cartan) generators 
- i runs from 1 to N — 1. Here Ejj is 2N by 2N matrix with one non zero element equal to 
1 and located at the position (i, j). The motivation for this choice is clear: the first 3 types of 
generators with opposite choice of sign between two E belong to Sp{N ) since are of the form 

a b 
c —a 

required for Lie algebra of Sp(N ) with the symplectic form choice we made, where a, b, c- are N 
by N matrices and b and c - symmetric. So, opposite choice of signs means that these generators 
are in orthogonal completion, i.e. in the coset algebra su(2N)/sp(N). Condition of unitarity for 
general g = exp (i n I (x)Ei^ now reads as n{* = nk, nfj* = rijj 1 , riff 1 * = rijf. It is clear that 
n 11 ,n I11 can be taken antisymmetric. 

Generators of sp(2P) / sp(P ) x sp(P) coset Lie algebra may be chosen in the form of 4P by 4P 
matrices, dividing them into 16 P by P block matrices.We define the action of the first invariant 
sp(P) subalgebra in the four blocks (1,1), (1, 3), (3,1), (3, 3) then the second sp(P ) acts in the 
blocks (2, 2), (2,4), (4, 2), (4,4). We require the generators of the coset Lie algebra just to be zero 
in these eight blocks and get the following basis: 

Ei,j+P -®7+3P;*+2P 
Ei+P„j Py+ 2P.i+3P 
Ei,j+3P + Pj+P,i+2P 
Ei+2P,j+P + Ej + 3P } i 

These are non zero in the remaining eight blocks, and are generators of sp{2P). The reality 
condition leads to the requirement nj = njp nni = niv- 

Calculation of bilinear terms may be done by expansion of exponent for coset group element 
up to the second order in fields and substitution of it into (2). For instance for this leads to 
Cq = —T7'(\ a \b)d tl n a d f j j rib 1 where A a , n a - a total set of Lie algebra generators and corresponding 
fields. Explicit calculation also shows that the term linear in h gives total derivatives, and may 
be omitted. In the same way, terms containing fields for Cartan generators n* drop out from the 
term proportional to h 2 in (2). So fields rii in both cases decouple and may be omitted. For the 
remaining fields we change notations and normalization: All case n/ —> n,nji —> ?n, CII case 
n/ —* n,nni —> m, and get the actions (3). 


Kj = 

ei; - 

pin _ 
piv _ 

E ij 


8.2 Projectors 


Here we present the projectors appearing in irreducible decomposition of tensor product of two 
antisymmetric representations. For projectors appearing in irrep decomposition of tensor prod¬ 
uct of two antisymmetric representations (highest weight ^ 2 ) for SU{2N) written in terms two 
antisymmetric tensors of rank 2 ajj and b^i, one can get using standard Yang tableau technique. 


P/M 

TV 

P/U+M3 


— ((ijjbki + aubjk T cikibij T djkbu ^ikbji djibik) 
-(2aijbki + 2akibij Y dkjbn -\~ ciikbji Y djib^ Y ctubkj') 

~{ a ijbkl CLklbij ) 


(20) 


9 



With more work one can get also projectors appearing in irrep decomposition of tensor prod¬ 
uct of two antisymmetric representations (highest weight ^ 2 ) for Sp(N ) written in terms two 
antisymmetric rank 2 tensors and bki traceless in the sense that J2 i,j a ij a ij = 0) where we 


chose for symplectic form <7 


0 1 

-1 0 


with four N x N block matrices. 




V, 


V1+V3 


^2,ii 

^2 

Vo 


g {cCijbkl T 0,1bjk ~h CLklbij “I - Cljkbil Oikbjl 0.jjbj}, 


( 21 ) 


+ 


N - 2 


4 / 

( a ^)[ifc] cr j7 + i a ^)[ii} a jk - (ab)[ji]<Jik + ( ab)[jk]au 
+ {ab)[ij](Tkl + { a b\kl}^ij'j ~k _ 2j ~ Vil&jk ~ &ij&kl) 


O'ijbkl o.kib,, T 


y 


N - 2 


{{ a b)(ii)Vjk + (ab) w au 


( a b)(ik) a ji i^b) (ji)&ik 


V 2 ^ .. {2 CLijbkl 2 CLklbij T 0>kjbil “k Ojkbj} T djlbik T Ojj bkj 


[(4 (ab)ki + 2 {ab)ik)(Jji - (4 (ab)u + 2 (ab)u)a jk 


N - 2 

+(4(a6) Zi + 2(ab) j i)cr i k - (4(o6) fcj + 2(ab) jk )cr i i - 4 (ab\ k i]crij 

(( ab )) 


+ 


fto’ikO’ji 2(Ju(jjk -|- 4(jjjO'ki) 


(N - 1)(N — 2) 

PikCTjl Bn&jk -|- BjiCik Bjk&U 

Cik&ji CuCTjk 4 “ Cji<Jik Cjkcru 

((ab)) 


N(N - 1) 


(CTikdjl CTj k (Tn ) 


where 


(a% 


/j 


y 


iV N 

'y ] O'ikbjicrki ; (( 0 ^)) = ^ ' (ab)ij(Tij, 
k,l= 1 i,j=l 


2(_/V — 0 ) + ( Q b)ji )) 


Cy = 


1 

2(iV - 2) 



(a% - y(( ab )) 
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